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Theory of quasiparticle transport in the mixed state of a d-wave superconductor is developed under 
the assumption of disordered vortex array. A novel universal regime is identified at fields above 
H* = c* H C 2(T/T C ) 2 , characterized by a, field-independent longitudinal thermal conductivity k% x . It 
is argued that this behavior is responsible for the high-field plateau in k xx experimentally observed 
in high-Tc cuprates. 



Properties of the mixed state of high-Tc cuprate su- 
perconductors keep surprising the scientific community 
in ways never imagined before. On one hand the collec- 
tive behavior of vortices produces a multitude of "vortex 
phases" attributable to an intricate interplay between the 
thermal fluctuations, dimensional crossover, and pinning 
forces on vortices. At low temperatures, on the other 
hand, unexpected and fascinating properties are observed 
related to the d-wave symmetry of the order parameter 
and the consequent relativistic "Dirac" spectrum of the 
low-energy quasiparticle excitations. 

A recent example of such an unexpected behav- 
ior is the experimental observation of the high-field 
plateau in the longitudinal thermal conductivity k xx 
in Bi 2 Sr 2 CaCu 2 8+5 (BiSCCO) and YBa 2 Cu 3 6+a; 
(YBCO) by Krishana, Ong and co-workers (TJQ] and by 
others |]-[^] ■ In particular it is found that after an initial 
steep drop at low fields (but well into the mixed state) 
k xx becomes field-independent above H* = c*H C 2(T/T c ) 2 
and remains so up to the highest attainable fields ~ 14T. 
The implication is that both electronic and phononic con- 
tributions to k xx separately become independent of field, 
with the initial drop attributed solely to electronic k xx 
^,|). These observations stand in a sharp contrast to 
the behavior of k xx (H) found in conventional s-wave su- 
perconductors H , where vortices are strong scatterers of 
both electrons and phonons. Understanding the physics 
of the plateau phenomenon therefore presents a consider- 
able challenge to the theory of quasiparticle transport in 
cuprates. Its broader significance lies in the fact that si- 
multaneous measurement of the field-dependent longitu- 
dinal (k xx ) and transverse (n xy ) thermal conductivities, 
now feasible experimentally, contains a great wealth of 
information on the quasiparticle dynamics and in princi- 
ple affords deduction of inelastic scattering rate, thought 
to be important for the pairing mechanism in cuprates. 

The initial interpretation of the plateau involved a 
field-induced transition to a fully gapped state, such 
as d x 2_ y 2 + id xy (or "d + id' v ), which would effectively 
freeze out the quasiparticle transport at low energies jjj . 
Laughlin presented a compelling theoretical argument 
in support of such a scenario with a weakly first or- 
der transition, but pointed out that opening of a gap 
with physically reasonable magnitude could not in itself 
account for the complete suppression of the electronic 



transport in BiSCCO at temperatures up to 30K. Fur- 
thermore, we find that sudden opening of a sizeable gap 
leads to a jump in k xx (H), as opposed to a kink observed 
in [0J|]. Experimentally, there appears to be little addi- 
tional evidence for the d + id' state, except perhaps for 
the apparent bound states found by scanning tunneling 
microscopy (STM) in the vortex cores of YBCO || which 
should not exist in a pure d x 2_ y 2 state [B. 

The existing theoretical treatments ]f^,|ll|-^13| so far 
side-stepped the important issue of the effect of mag- 
netic field on the quasiparticle mean free path (MFP), 
determined primarily by scattering of quasiparticles by 
the Abrikosov vortices. Detailed theory of the MFP due 
to vortices in s-wave superconductors was developed long 
ago by Cleary jl(|. The ci-wave problem, however, is fun- 
damentally different [9] and, except for heuristic treat- 
ments given in Refs. ]14|Jl^ |, no analogous theory exists 
at present. Here we develop such a theory and show 
that field-independent longitudinal thermal conductiv- 
ity arises quite naturally in a pure d x 2_ y 2 state above 
a crossover field H*(T). This behavior reflects an ex- 
act compensation between the enhancement of the quasi- 
particle density of states (DOS) in the presence of non- 
uniform superflow (the "Volovik effect" ) |l6| and the con- 
comitant reduction in the quasiparticle mean free path i 
due to increased scattering from vortices in a disordered 
vortex array. The limiting high- field value of k xx (H) is 
universal in the similar sense as the impurity induced 
universal microwave conductivity <t(lu — > 0) predicted by 
Lee [jlTj . The approach to this limiting value depends on 
the distribution of vortices in the sample and will there- 
fore be material and sample dependent. 

The basic physics of the new universal regime can 
be understood from the following simple argument. In 
the elementary Boltzmann-type treatment the electronic 
thermal conductivity of a normal metal is given by 

K £ xx = ^V F C V £, (1) 

where c v is the electronic specific heat and £ is the MFP. 
We now argue that this relation remains valid in the 
mixed state of a d-wave superconductor at fields above 
H*. As first pointed out by Volovik |l6| ], the superfluid 
velocity field around vortices induces a Doppler shift in 
the excitation spectrum of quasiparticles, which in turn 
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leads to a finite DOS at the Fermi level. For H > H* 
the system behaves effectively like a normal metal and 
we expect ([j]) to hold. The residual DOS scales as 
Nh(0) ~ \H and gives rise to the well known low tem- 
perature specific heat mM 



koN F Ty/H/H c: 



(2) 



where Np is the normal-state DOS, and kg is a constant 
of order unity. This type of behavior is indeed observed 
experimentally fig] ]. 

We now estimate the vortex contribution in to the 



total MFP t 



The central assumption 



we make is that in the regime of experimental interest 
the vortex lattice is disordered in the sense that it pos- 
sesses no long range translational order, and will there- 
fore scatter quasiparticles (fl|. At low energies £h is 
dominated by the quasiparticle scattering from the su- 
perfmid velocity field v s (r). Scattering from the vortex 
cores is down by a factor (£/a„) 2 = H/H c2 <C 1 p2[ , 
where a v — £ yj H C 2 / H is the average inter- vortex dis- 
tance, £ — vf/ttA is the coherence length, and A is the 
maximum gap. In the (disordered) vortex lattice v s (r) 
varies on the length scale set by a v . Since a v is the only 
relevant length scale in the problem |23] , one expects on 
general grounds that 
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where k\ is a field-independent constant. This conclusion 
is indeed confirmed by an explicit calculation of the quasi- 
particle propagator outlined below, as well as by a calcu- 
lation of the vortex transport scattering cross-section otr 
p2f , carried out along the lines of the classical treatment 
by Cleary [io| ]. We note that Eq. (||) is also consistent 
with the result of Ref. jl5| based on a heuristic argument 
involving the Andreev reflection. 

In sufficiently strong fields Eq. (|J) implies that In "C £o 
pij and the total MFP will be dominated by vortices: 
£ « in- On substituting ^ and (^) into (|l|) we arrive at 
the desired result that for H ^> H* , k xx will approach 
the field-independent universal value 



K e x */T = k'irN F v F /3A , 



(4) 



with k' = kok±. Evaluation of k 1 below shows that K e Jf is 



In the absence of field the quasiparticle propagator is 
a 2 x 2 matrix in the Nambu space (taking h = 1) 



Go(w,k) = 



r 3 ek + riA k 



(•5) 



where fj are the Pauli matrices. In intermediate mag- 
netic fields the coupling of quasiparticles to the superflow 
around the vortices is well described by the semiclassical 
replacement oj — > lo — k • v s (r) in Eq. (||) |TJ,|2^]. In 
a simple London model, which will be sufficient for our 
purposes, the superfluid velocity field is given by p8| 
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Here A is the London penetration depth, z is a unit vector 
along the field direction, and {Ri} denotes vortex posi- 
tions. In a disordered vortex array, on the length scales 
large compared to a v , propagation of quasiparticles will 
be described by the Greens function averaged over {R;}: 



G( W ,k)HG ( w -k.v s ,k)) {Ri} . 
If we now define a probability density 

V(r ] ) = (S[r 1 -k.v s (r)}) {Ri} , 

we may rewrite (Q) as 



G(w,k) = / dr)V(r))G (LJ-r),-k), 



(J) 



(8) 



(9) 



and all the information on the vortex array is now en- 
coded in V{r\). Making the most natural ansatz that the 
vortex positions are random and uncorrelated, it can be 
shown p| that V(r)) = (27rcr 2 f )- 1 / 2 e-'' 2 / 2<J « with 



^ H = k a k f3 (v:(r)v^r)) {Ri} ^^(^- 



Aq hiK. 
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Here k — X/£, and the last equality follows from Eq. (||). 

Evaluation of the propagator G is somewhat compli- 
cated by the fact that the 77-integral in Eq. ([)]) cannot 
be expressed in terms of elementary functions. This dif- 
ficulty can be avoided by taking V to be a Lorentzian 
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ay) rather than a Gaussian, 



identical to the "universal" thermal conductivity Knn due u - 1 ,1 • , ,• • , , j u , • 

. rfrSl FTlH j which case the integration is elementary and one obtains 



to impurities predicted by various authors [|12|,g5| and 
observed experimentally by Taillefer and co-workers [^6| . 

The argument based on Eq. (|l|) captures the essential 
physics of the new universal regime and provides a very 
natural explanation of the plateau phenomenon observed 
in cuprates jj]-^ in terms of fundamental properties of 
the Dirac fermions. We now present a more rigorous 
treatment of k% x (H) , based on the Kubo formula for the 
heat current response. Besides supplying the unknown 
constant fc', such calculation provides the necessary con- 
fidence in our result (Q) and helps understanding the ap- 
proach to the universal limit with the increasing H . 



G(u>, k) = Gq(uj — i<jH,k). 



(11) 



Scattering from vortices therefore results in a self-energy 
correction to the Greens function, which, for the partic- 
ular case of Lorentzian distribution, is simply a constant 
S^f(w, k) = i<7H- Had we kept the Gaussian distribu- 
tion, the self energy would be similar at low energies, 
but would tend to zero for \uj 2 — e\ — 3> o 2 H . This 
additional structure results in slight mathematical com- 
plications but does not change the result for k% x (H) qual- 
itatively. We therefore choose to proceed with the simpler 
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form (|11|) and treat it either as a reasonable low energy 
approximation to Eq. @ or as a slightly different physi- 
cal problem with a vortex distribution {R^} resulting in 
the Lorentzian distribution V(rj). 

Eq. (|ll]) may be used to calculate the field-induced 
DOS, N h (lu) = -(2vr)- 1 X; k Tri(w,k), where i(w,k) = 
2ImG(o;,k) is the spectral function. One obtains 
Nh(0) oc an oc \/H, in agreement with Volovik's result 
fliql . Interpretation of an as the scattering rate due to 
vortices allows estimation of the vortex MFP £jj = v/cth, 
where v is the average quasiparticle velocity. Taking 
v = ,/v~fva fll|] we find £jj = o, v (8va/vf In k) 1 / 2 ~ 0.5a„ 
(for Vf/va = 7 and k = 70), in agreement with our naive 
estimate (||). We finally note that since the scattering 
rate an is proportional to ~ y/nZ (with n v = H/&q the 
vortex density), Eq. (11) will be in general difficult to 
derive from conventional diagrammatic techniques which 
typically employ low-density approximation and there- 
fore yield scattering rates proportional to n v . Physically, 
the yfnZ non-analyticity reflects the long ranged 1 jr na- 
ture of the v s field associated with a single vortex. 

We are now in the position to evaluate the thermal con- 
ductivity tensor from the standard linear response the- 
ory. Following Ambegaokar and Tewordt |29| , we find for 
a d-wave superconductor in 2D, 
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where (neglecting vertex corrections) 



d 2 k 

(27T) 2 



fc l ^Tr{f 3 i(w,k)f 3 i(w,k)}. (13) 



At temperatures low compared to Ao it is permissible to 
linearize the quasiparticle spectrum in the Dirac cones 
near the four nodes of the gap function Ak = Ao cos(2#) 
0,0. Taking £k = v F ki and A k = in the new 

coordinate system (fci,fc2) with the origin at the node, 
one can explicitly perform the integral in ( |l3| ) to obtain 



K xx (lo) 
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Here va — 2Ao/fcF and we have replaced an by the to- 
tal scattering rate a = <jq + ajj, with <7q describing both 
elastic and inelastic processes in the superconductor at 
H = 0. The proper description of these processes would 
presumably require a self energy £o(u>, k) which depends 
strongly on both of its arguments. However, absent the 
detailed microscopic theory of such processes, we follow 
fjof and simply model their effect by a phenomenologi- 
cal constant ao, which we interpret as a process-specific 
average of So (a;, k). 

For arbitrary a and T the expression (12) for k% x must 
be evaluated numerically. However, the leading terms 
may be readily obtained in a perturbative expansion. We 
obtain, for a ^> 2T (high field regime): 
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and for a <C 2T (low field regime): 



v F N F 
2A n 



9C(3)T ln2cr 
4 a + ~T 



(15) 
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The first term in Eq. ( |T5| ) reproduces Eq. (||) for k' — 
1/2 and represents the universal thermal conductivity 



iQ |Tjj25[], which does not depend on the de- 
tails of the vortex distribution, as long as there is no 
long range order. This resultjustifies the usage of the 
simple Boltzmann approach (|lj) in the mixed state of a 
d-wave superconductor. The leading correction in ( |l5| ) 
is non-universal; e.g. one may show that the power will 
change to (T/ct) 4 for a Gaussian V '. 

Eqs. (|l5|) and ([if]) provide a simple tool for extracting 
the zero-field scattering rate ao from the experimental 
data. Assuming that the phonon contribution ft?Lto the 
thermal conductivity field-independent (M, the 

total drop in n xx {H) between H — and the plateau is 
simply related to ctq by 



a = 



ciT(l + Sk xx ) 
c 2 T(fc ra )- 1 / 2 , 



do «T 
to »T 



(17) 



Here 5k xx = [k xx (0) - k xx (oo)]/ Kq , k xx (oo) denotes 
the plateau value, c\ — 27C(3)/47r ~ 2.58 and c 2 = 
y/7w 2 /15 ~ 2.15. Applying these relations to the data on 
underdoped YBCO § and using n e m /T = 0.019W/K 2 m 
p6]| , we obtain ao ~ 0.4meV at 10K, increasing to about 
7meV at 50K and 19meV at 60K (just below T c = 63K). 

To illustrate the behavior of k% x (H) over the wide 
range of fields and temperatures in Fig. |l| we have eval- 
uated Eq. ([l2]) with kernel ([l4|) numerically, adopting 
a phenomcnological expression for the scattering rate 
(To/Ao = 70+73(17 A ) 3 . Here 70 represents the residual 
scattering rate due to impurities and the 73 term models 
the inelastic scattering rate |^] which is known to col- 
lapse rapidly below T c [Q. Similarity of Fig. |l| to the 
experimental data on the underdoped YBCO M is strik- 
ing: the positive curvature of n xx (H) is enhanced with 
decreasing T and the curves approach a field-independent 
value at high fields. The characteristic crossing of the 
curves near T/Ao ~ 0.15, which reflects the crossover in 
ao(T) from inelastic scattering at high T (solid lines) to 
elastic scattering at low T (dashed lines) is also present. 

In some samples of BiSCCO, the plateau is reached 
abruptly, with a discontinuity in the slope of n xx (H) at 
H* |l],p|)[|) indicative of a phase transition. This "kink" 
behavior is not reproduced by the present simple model, 
which, however, can be presumably appropriately modi- 
fied once the nature of the phase transition is understood. 
We take a position here that it is the plateau phenomenon 
which is universal and significant (since it is observed in 
more than one compound) while the approach to it is a 
material-specific issue of secondary importance. 
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FIG. 1. Normalized electronic thermal conductivity 
n xx (H) as a function of H at selected temperatures as cal- 
culated from Eq. (p^). Parameters used are H C 2 = 160T, 
k = 70, 70 = 0.1 and 73 = 4. 



Aubin and co-workers |3| further report that the be- 
havior of k xx (H) in BiSCCO qualitatively changes be- 
low IK, in that it increases with magnetic field, approxi- 
mately as k xx (H) ~ y/H. From Eq. (|]) the most natural 
interpretation is that, below IK, £ becomes independent 
of H . There are two possible reasons for this. The vor- 
tex array may order at low T, in which case the Bloch 
theorem prevents quasiparticles from being scattered by 
vortices Uj . The second possibility stems from the com- 
putation of the transport scattering cross-section of a d- 
wave vortex 22 1, which indicates that Otr(w) vanishes as 



for quasiparticle energies u> < Ao/(A;.fA). Thus, at 
low energies, vortices become "transparent" , and even 
disordered vortex array will not scatter quasiparticles. 

An immediate consequence of our picture is that all 
of the high-T c compounds with the d x 2_ y 2 gap should 
exhibit the plateau phenomenon with the universal value 
K xx /T — ttVfNf/6Ao in sufficiently high fields. We pre- 
dict that the existence of plateau will correlate with the 
absence of long range order in the vortex array, and will 
fade away at the lowest temperatures. The formalism 
developed here also permits calculation of the Hall con- 
ductivity n xy . However, this constitutes a more compli- 
cated problem since one has to consider the particle-hole 
asymmetry in the vortex scattering rate ]lC| ] as well as 
corrections to the linearized Dirac spectrum jnj. 
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Ando, A. V. Balatsky, K. Behnia, A. J. Berlinsky, P. 
J. Hirschfeld, C. Kallin, A. J. Millis, L. Taillefer and Z. 
Tesanovic for valuable discussions. This work was sup- 
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at the Aspen Center for Physics. 
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